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ABSTRACT 

We present new models of the deuterium chemistry in protoplanetary disks, 
including, for the first time, multiply deuterated species. We use these models 
to explore whether observations in combination with models can give us clues 
as to which desorption processes occur in disks. We find, in common with other 
authors, that photodesorption can allow strongly bound molecules such as HDO 
to exist in the gas phase in a layer above the midplane. Models including this 
process give the best agreement with the observations. In the midplane, cosmic 
ray heating can desorb weakly bound molecules such as CO and N 2 . We find 
the observations suggest that N 2 is gaseous in this region, but that CO must be 
retained on the grains to account for the observed DCO + /HCO + . This could 
be achieved by CO having a higher binding energy than N 2 (as may be the case 
when these molecules are accreted onto water ice) or by a smaller cosmic ray 
desorption rate for CO than assumed here, as suggested by recent theoretical 
work. 

For gaseous molecules the calculated deuteration can be greatly changed by 
chemical processing in the disk from the input molecular cloud values. On the 
grains singly deuterated species tend to retain the D /H ratio set in the molecular 
cloud, whereas multiply deuterated species are more affected by the disk chem- 
istry. Consequently the D/H ratios observed in comets may be partly set in the 
parent cloud and partly in the disk, depending on the molecule. 

Subject headings: circumstellar matter — ISM: abundances — ISM: molecules 
- solar system: formation — stars: formation — stars: pre-main-sequence 
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Introduction 



Observations of deuterated molecules have long been used to trace the physics and 
chemistry of interstellar clouds. The elemental abundance of deuterium relative to hydrogen 
is ~ 10~ 5 , but the relative abundances of deuterated to non-deuterated molecules can be 
much higher, especially in cold regions. Deuterium is transferred from ions such as H 2 D + 
and CH2D + via gas phase ion-molecule reactions. The reverse reactions, removing the deu- 
terium, are inhibited at low temperatures, due to the difference in the zero point energy 
of the deuterated species to its non-deuterated equivalent. Another route to high levels of 
deuteration is via reactions between species accreted on to the surfaces of dust grains. These 
proce sses can lead to molecular D/H ratios of > 0.01 in cold molecular clouds (e.g. iTurner 
20011 ) . High levels of deuteration are also seen in hot cores, protostellar disks and comets. 



Recently, multiply deuterated molecules have been observed for the first time with high 
levels of deuteration, e.g. NHD2/NH3 ra nges from ~ 0.005 (in the cold core, L 134N) to ~ 
0.03 (in the low mass protostar 16293E) faoueff et allbood; lLoinard et al.llioOlh ND-J NfL 



Lis et al. I2OO2I ; van der 



~ 10^ 3 (in NGC133 IRAS4A and Barnard 1 
D2CO/H2CO ~ 0.01 - 0.4 (in cold cores and low mass protostars 



ak et al.ll2002h. and 



Loinard et al. 



2002 



Bacmann et al.l 120031 ) . High abun dances of deuterat e d methanol and formaldehyde have 



been seen in star forming regions (IParise et al.l |2002| ; lLoinard et al.l |2002| ; iBacmann et al. 



20031 ) and are attributed to the thermal desorption of grain mantles accumulated during an 
earlier low temperature phase, demonstrating the importance of grain surface reactions in de- 
termining deuteration. In order to achieve such high molecular deuteration ratios, very high 
atomic D abund ances are required during the period w hen ice mantles can ac cumulate, with 
D/H ~ 0.2 - 0.3 JCaselli et al.ll2002l : IParise et al]l2002h . The recent models of iRoberts et al. 



( 120031 ) have shown that the inclusion of HD^ and D3" can drive up the atomic D /H ratio to 
a level where grain surface reactions can account for the observations of multiply deuterated 
H 2 CO and CH3OH. 



A. few deuterated molecules have now been observed in protostellar disks. Ivan Dishoeck et al 

J200J observed DCO+ in TW Hya with DCO+/H CO+ = 0.035. This mol e cule was also ob- 
served in DM Tau, with DCO+/HCO+ = 4 x 10" 3 fiGuilloteau et alJbood ). ICeccarelli et~aL 



( 120041 ) observed H2D + in TW Hya and DM Tau and found fractional abundances of 4 x 10~ 10 
and 7 xlO -10 respectively. They argued that the emission arises in the midplane, where few 
other molecules exist because of efficient freezeout, and th us observations o f ^D" 1 " provide a 
means of determining the ionization fraction in this region. ICeccarelli et al.l (120051 ) presented 
observations of HDO in DM Tau with a fractional abundance of ~ 3 x 10~ 9 . The emission 
arises in a region with temperatures below 25 K, where HDO is expected to be frozen onto 
dust grains. The presence of HDO emission shows that an efficient non-thermal desorption 
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process must be acting, which iDominik et al.l (120051 ) suggest to be photodesorption. The 
detection of HDO and H 2 D + in DM Tau has been disputed by iGuilloteau et al.l (120061 ). 



The similarity between the molecular D/H ratios observed in comets and those seen 
in molecular clouds is one reason why it has b een su gges ted that comets may be com - 
posed of interst e llar m aterial, e.g. iMeier et al.l (ll998aU bl); iBockelee-Morvan et al.l (119981 ); 
Eberhardt et al.l (119951 ). This would mean that the molecules had experienced little or no 
chemical processing in th e protosolar disk, y e t mod e ls of deuterium chem istry in disks, e.g. 
Aikawa fc Herbstl (Il999ah : lAikawa fc Herbstl (1200 lh : lAikawa et al.l (|2002|) show that exten- 
sive chemical processing can occur, resulting in D/H ratios that are very different from the 
ones set in the parent molecular cloud. These disk models have considered the chemistry of 
singly deuterated species only and did not include grain chemistry, there by ignoring a po 



tentially important process for determining the D/H ratios. The model of I Aikawa fc Herbst 



( 1l999al ) considered the evolution of the chemistry in a cold cloud core, through the infall 
stage and into the protostellar disk. They found, in the region of the disk where comets are 
expected to form, that the molecular deuteration is affected by chemical processing in the 
disk. Therefo re the abundances of mo lecules in comets does not just reflect interstellar cloud 
abundances. lAikawa fc Herbstl (120011 ) considered only the disk chemistry, b ut used a disk 



model that was isothermal at any given radius. This was improved upon by lAikawa et al. 



( 120021 ) w ho used temperature a n d density distribu t ions derived from the hydrodynamical 
models of Id'Alessio et al.l (11999I ). lAikawa fc Herbstl (120011 ) also ignored grain reactions, but 
they found good agreement with the observed DCN/HCN and DCO + /HCO + ratios (~ 0.01). 

Models have shown that in highly depleted regions, such as the d isk midplane, the 
deute ration of Hjj~ is efficient and Djj~ can become the dominant isotopomer (jCeccarelli fc Dominik 
20051 ). This has consequences for the d euteration of o t her molecul es since they form from 



reactions of H^~, H?D + , HD 



and Djj~. 



Roberts et al.l (120021 . 120031 ) demonstrated that the 



molecular D/H ratios in dense molecular clouds can be significantly enhanced by the inclu- 
sion of HD^ and D^ in the models. This effect is greatest in regions where molecules are 
heavily depleted by freezeout onto grains. Since disks are very dense and therefore experi- 
ence a high degree of depletion, especially in the midplane, it seems likely that HD^~ and Djj~ 
could have a significant effect on the chemistry. 

Given the importance of deuterium in tracing the origin of solar system bodies and for 
understanding the thermal history of star forming regions, together with the evidence that 
multiply deuterated forms of Hjj~ are important for determining the molecular deuteration 
and the suggestion that grain surface reactions may be important in determining the D/H 
ratios, it seems timely to take another look at the chemistry of deuterium in protoplanetary 
disks. Here we present a model of the chemistry of a T Tauri disk that includes multiply 
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deuterated species, combined with grain surface chemistry. We consider how different des- 
orption processes can affect the column densities and compare our results with observations 
with the aim of better understanding the processes that are at work in disks. We also use 
our models to consider how processing in disks may affect the molecular D /H ratios observed 
in comets. 



2. The models 

In this section we discuss the processes that are included in our models. We wish to 
investigate how different assumptions about the chemistry, in particular about desorption, 
can affect the calculated abundances. To do this we have chosen one particular physical 
model to provide the density and temperature distributions. We consider three desorption 
processes, thermal desorption, desorption due to cosmic ray heating of grains (CRH) and 
photodesorption. All models include thermal desorption. Model A excludes non-thermal 
desorption, Model B includes CRH, Model C includes CRH and photodesorption, and Model 
D includes photodesorption but not CRH. Our models assume that the disk is static, and 
that no mixing processes are acting. 



2.1. The chemical model 



Our basic chemical model is taken from the UMIST database RATE95 (Millar et al. 



1998). The deu t erate d net work is derived from this database using the techniques described 



in 



Millar et all (119891 ) and iRodgers & Millar! (119961 ). In the gas phase formation of mono- 
deuterated species the deuterium is transfer red by reaction of H 2 D + , CH 2 D + and C2HD~ 



The rates for these reactions are taken from iRoberts &: Millar! (120001 ) . For multiply deuter- 



ated species the reactions of H D+, Dj~, 



CHD+, CDt and C2HD + are important. Rates for 



these reactions are taken from IRoberts et all (120041 ) . Other reactions involving deuterated 
species are assumed to have the same rates as their non-deuterated counterparts. For reac- 
tions involving complex species, the branching ratios must be determined. In the absence 
of any information on the reaction mechanism, this is done on a statistical basis, with the 
exception of reactions where the same functional group appears as both a reactant and 
a product. In this case it is assumed that this group is preserved through the reaction 



(IRodgers &: Millarlll996l ). For example, in the dissociative recombination of CH2DOH + only 
CH 2 DOH is formed and not CH 3 OD. 



Our reaction network links 227 gas phase species (of which 115 are deuterated) and 
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91 grain species (44 deuterated) by 9489 reactions. In order to reduce the time taken for 
the model to run, not every possible multiply deuterated species is included e.g. multiply 
deuterated forms of methanol are not included. 



2.2. The recombination of N 2 H 



In RATE95 it is assumed th at N 2 H + recombines with electrons to form N 2 and H only. 
However, recent laboratory work (IGeppert et al.ll2004l ) has shown that this is in fact a minor 
reaction pathway and instead 



N 2 H+ + e 



NH + N 



is the main route. This result has consequences for the chemistry. iRoberts et al.l (120031 ) 
found that the new branching ratios do not affect the calculated abundances in gas phase only 
models, but do have significant effects in models that include freezeout. At low temperatures 
when freezeout is included, the destruction of N 2 H + by reaction with CO becomes less 
significant as the CO is depleted, and so Reaction [T] gains in importance. The formation 
of NH from reaction [T] results in the removal of nitrogen from the gas in cold regions, since 
any NH that accretes onto the grains will be rapidly hydrogenated to form NH3, a molecule 
with a high binding energy that is therefore not easily desorbed, unlike N 2 which is very 
volatile. Consequently, the new reaction pathways can reduce the abundance of nitrogen- 
bearing molecules in the gas. This reduction is less significant in regions where the grain 
temperature is higher than about 20 K, where NH can be thermally desorbed before reacting, 
and hence the formation of NH3 ice is less efficient. 



In disks, ICeccarelli fc Dominiki (120051 ) showed that N 2 is an important means of control- 
ling the degree of deuteration. The presence of N 2 will prevent the transfer of deuteration 
along the chain of isotopes from H3 to D3 and reduce the abundance of D3 by destroying 
the less deuterated isotopomers before they have a chance to form Djj". The new reaction 
pathway means that the abundance of N 2 will be reduced in the midplane, and hence the 
degree of deuteration of H^ will be higher. 



2.3. Freezeout and desorption processes 

Given the high densities and cold temperatures found in much of the outer disk, molecules 
that hit a grain are likely to stick to it efficiently. We assume that all species freezeout at 
the same rate with a temperature independent sticking coefficient of 0.3. For positive ions 
the freezeout rate is increased slightly, since grains are likely to have a negative charge 
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( jUmebayashi fc Nakanol Il980l ) and therefore there could be a stronger attraction between 
the positive ions and the negatively charged grains. The freezeout rate of positive ions is 
assumed to increase by a fact or C = 1 + 16.71 x 10~ 4 /(aT g r ) where a is the grain radius and 



T gr is the grain temperature (jUmebayashi &: Nakanol Il980l ). Ions are assumed to recombine 



on the grain surfaces in the same way that they do when reacting with electrons in the gas, 
e.g. in the gas phase 

HCO+ + e — > H + CO (2) 



and on the grain 



HCO 



gr 



H : gr + CO : gr 



(3) 



where H: gr and CO: gr represent H and C O on the grains. The exceptions are N2H + and 
N2D" 1 " where we follow iRoberts et al.l (12004 ) in assuming that these ions form N2 rather than 
NH and ND when they freezeout. All species are assumed to freezeout, except for He which 
has a very low binding energy and is therefore easily thermally desorbed even at very low 
temperatures. Any He + that hits a grain is neutralized and returned immediately to the gas. 

We include thermal desorption and desorption due to cosmic ray heating of grains. 
Cosmic ray heating is able to maintain a low level of some volatile molecules such as CO and 
N2 in the cold midplane of our models. These molecules can destroy ions such as Hjj" and its 
deuterated isotopes, and produce ions such as HCO + and DCO + in the midplane and thereby 
affect the ionization level in this region. We have therefore run models without cosmic ray 
heating to determine how the inclusion of this process will affect molecular abundances and 
the i onization level in cold regio ns of the disk. The rates for cosmic ray heating are taken 
from lHasegawa fc Herbstl (119931 ). using updat ed binding energie s for some species, notably 
CO , for which we use the value determined by lOberg et al.l (120051 ) . Table [1] gives the binding 
energies (Ed) used in our models. 



Thermal desorption rates are calculated from 



k 



therm 



u e 



-E D /T g . 



(4) 



where T gr is the grain temperature, and vq is the frequency of oscillation between the ab- 
sorbate and the surface given by 



V(2n s E D / 



7r 2 m) 



(5) 



where n s is the surface density of sites (~ 1.5 x 10 15 cm 2 ) and m is the mass of the accreting 
species. 
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2.4. Photodesorption 



Ceccarelli et al.l (120051 ) recently observed a high column density of HDO in the disk 
around DM Tau. They found N(HDO) ~ 1.6 x 10 13 cm -2 in the outer disk, where the 
density is ~ 10 6 cm -3 and the temperature is < 25 K. The emission comes from above 
the midplane and corresponds to a relatively high fractional abundance of ~ 3 x 10~ 9 . At 
these densities and temperatures, water is expected to be completely removed from the gas 
by accretion onto grains, where it will remain in the absence of a non-thermal desorption 
process. Cosmic ray heat i ng is n ot efficient enough to remove such a strongly bound molecule 



as water. IDominik et al.l (120051 ) suggested that photodesorption, arising from the action of 
the interstellar far-UV field could be efficient enough to retain water vapor in the gas phase 
in the cold outer disk, and so account for the observed column densities. They found that 
photodesorption can maintain a layer of water vapor above the midplane with a fractional 
abundance of ~ 3 x 10~ 7 . The average calculated value of N(H20) in their model was ~1.6 
x 10 15 cm -2 , and using this combined with the observed N(HDO) they deduce HDO/H 2 
= 0.01 in DM Tau. 

We have previously investigated the effects of photodesorption on the chemistry of disks 
( jWillacy fc Langerll2000l ) and found that it can retain high abundances of many molecules in 
a layer above the midplane, even in cold disks. Here we revisit the idea of photodesorption 
and investigate not only how it affects the abundances but also how it impacts the deuteration 
of molecules (Models C and D) 



The rate of photodesorption is given by 



k, 



pd 



GnYe- 1 '^ 



ira 2 n n 



(6) 



where Go is the radiation field in units of the Habing field (10 8 photons cm _2 s _1 ), Ay is 
the visual extinction and Y is the photodesorptio n yield. We use the temperature depen- 
dent value of Y as determined experimentally by IWestley et al.l (119951 ); From Figure 2 of 
Westley et al.l Y = 0.003 molecules per photon at T < 50 K and rises to 0.075 at T = 100 
K. 

We include the effects of desorption caused by photons from the i nterstellar radiation 
field, the stellar field and from the cosmic ray induced photon field ( jPrasad fc Tarafdar 
19M). 
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2.5. Grain surface chemistry 



Previous models of deuterium chemistry in disks have ignored grain surface reactions 
and thereby have excluded a potentially very important contribution to the molecular D/H 
ratios. The inclusion of grain chemistry is problematic, as the simplest way is to use the rate 
equation metho d, but this g i ves ver y different surface abundances than the more e xact Monte 



Carlo model. ICaselli et al.l (120021 ) found that using the H-atom scan rates of iKatz et al. 



(119991 ). which are much lower than usually assumed, produces grain mantle abundances in 
better agreement with Monte Carlo models than the results from standard ra t e equ ation 
models. The work of Katz ETA al. has been disputed, e.g. by iHornekaer et al.l (120031 ) who 
find that H 2 forms very efficiently on grains, possibly by H-atoms tunneling through barriers 
on the grain surface. However, because of the advantage in going some way to correcting 
the short-comings of the rate equation method we use the Katz wt al. scan rates here, and 
assume that the scan rate for D-atoms is also slow. 

We assume that reactions on grain surfaces can occur only if one of the reactants is 
an atom. All other specie s are assumed to be immobile. The reaction set is taken from 
Hasegawa fc Herbstl (119931 ) with the addition of the equival ent deuteration reac tions. Acti- 
vation barriers for the deuterium reactions are taken from ICaselli et al.l (120021 ) . The rates 
are temperature dependent and the reaction rate between two species, i and j, is given by 



k ij = ««(*, 7 +tj 1 )/(N s n d ) 



(7) 



where = e~ Ea ' kTgr (E a is the activation barrier to the reaction), N s is the total number 
of sites on the surface of the grain and is the grain number density. The rate at which % 
and j scan the surface of the grain is given by t% and tf 



t: 



z/ e 



(-0.3E D /T gr ) 



(8) 



where Ed is the binding energy (from Table [T]) and vq is given by Equation [5j 



Important parameters are the binding energies of H and D atoms and their reaction 
rates on the grains. For the formation rate of H 2 we have used the work of ICazaux &: Tielens 
(120021 . 12OOJ) who developed a model of this process that fits the experimental data for the 
reaction of two hydrogen atoms on silicate and amorphous carbon grains. They find that H2 
formation can be efficient at temperatures up to 500 K on these surfaces. This model takes 
into account the possibility of both chemi- and physi-sorption and assumes a high binding 
energy of 600 K for atomic hydrogen atoms adsorbed onto a silicat e surface. This value of 
Ed(H) is rather higher than the 350 K lTielens fc Allamandolal (119871 ) which we have assumed 
in our previous work, and which results in a greatly reduced H2 formation rate at T > 15 K. 
For deuterium atoms we follow ICaselli et al.l J2OO2I ) and take E D (D) = E D (H) + 21 K (21 



- 9- 



K is the zero point energy difference between the hydrogen and deuterium atoms). Recent 
calculations of the bind ing; energy of H-atom s adsorbed o nto water ice have found values 
ranging from ~ 400 K (Ul-Halabi et ajjooj ) to - 575 K JPerets et alil2005h . Our choice 
of 600 K is therefore a little high, but consistent with the larger of these calculated values. 
Given the uncertainty in Ed(H) we compare the results for Ed(H) = 600 K with those for 
£'d(H) = 350 K in Section f3. 1.51 to determine the effects of our choice on the chemistry. 



2.6. Ionization processes 

Ionization in the disk can arise from several sources: 



1. UV photons 

In the surface layers ionization by UV photons from both the stellar and interstellar 
radiation fields is important. We assume that the stellar photons travel out radially 
and that the interstellar photons hit the disk vertically. The two fields are combined to 
give an overall UV flux at each position in the disk. The strength of the stellar UV field 
for a T Tauri star has been estima t ed as 10 4 times the interstellar radiat ion field (Go) 



at 100 AU by iHerbig fc Goodrich! (119961 ). More recent observations by iBergin et al. 
( 120031 ) have estimated a much lower value of a few hundred times Go- We have chosen 
a value of 500 Go & t 100 AU. In common with lAikawa et al.l (120021 ) we assume that the 
stellar field does not dissociate CO and H 2 , although th ese molecu l es are dissociated 
by the interstellar radiation field. We use t he approac h of iLee et al.l (119961 ) to describe 
the self-shielding of both molecules. The ILee et al.l model is based on a slab model 
of a molecular cloud and provides data for the shielding due to H 2 , CO and dust 
as a function of column density. The linewidths assumed are 3 kms -1 - much larger 
than observed in disks, whe re the velocity dispersion in the outer disk is almost thermal 



Guilloteau fc Dutre 



vl il998h . 



To take account of this we have followed Aikawa fc Herbst 



1999 



b) in scaling the column densities in Table 10 of Lee et al. by c s /3 kms 1 where 



is the sound velocity. This scaling factor is only required for H2 since CO dissociation 
lines are broader due to predissociation. 

HD does not self-shield. However some of its lines do overlap with those of H 2 al - 
lowing shielding by the H 2 lines to reduce its photodissociation rate. iBarsuhnl (119771 ) 
estimated that the HD photodissociation rate will be reduced by 1/3, assuming that 
the overlapped HD lines are totally shielded by H 2 . We adopt this factor and scale it 
where necessary to take account of the region where H 2 does not fully self-shield. 



2. Cosmic rays 
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Cosmic rays can produce ionization in the disk if the surface density is low enough 
for them to penetrate i.e. less that 150 g cm -2 (jUmebayashi Nakandll98ll ). This 
is always the case in the region of the disk considered here. Cosmic rays can cause 
ionization both directly and indirectly by producing photons from their interaction 
with H2. The rates of these processes are taken from the UMIST ratefile with an 
assumed cosmic ray ionization rate of 1.3 x 10~ 17 s _1 . 

3. Radioactive isotope decay 

The decay of radioactive nucleides is an additional source of ionization. 26 Al can decay 
to form excited 26 Mg, which in turn decays by either positron decay or by electron 
capture. We include ionization due to these processes with a rate: 



Caz = 6.1x10- 18 s- 



(9) 



( lUmebayashi fc Nakanolll98ll ) 



4. X-rays 



T Tauri stars are strong X-ray emitters. Ilgea fc Glassgoldl (Il999l ) modeled the effects of 
X-ray ionization in a disk and found that it can be effective in the surfac e layers where 
the attenuation length of X-rays is very small. lAikawa fc Herbstl (1200 ll ) included the 
effects of X-rays in their model of the chemistry in the outer disk, but found that they 
did not significantly affect the results for this region. Based on this we have chosen to 
ignore the effects of X-rays in these models. 



2.7. Input abundances for the disk model 

We assume that the material which is incorporated into the disk has first been processed 
to some extent in the parent molecular cloud. We therefore use the abundances produced 
by a molecular cloud model that has been run for 1 Myrs as inputs to the disk model. 
The elemental abundances used as inputs to the cloud model are listed in Table EJ We 
assume that initially all hydrogen is molecular and that the deuterium is in HD. All other 
elements are in their atomic form, with the exception of carbon which is ionic. The molecular 
cloud model is run at a density of 2 x 10 4 cm -3 , T = 10 K and a visual extinction of 10 
magnitudes. It includes all the processes that are included in the disk model, with the 
exception of photodesorption and ionization by the decay of radioactive nucleides. The 
molecular abundances which are used as inputs to the disk model are given in Table [3] For 
comparison we also include the molecular abundances observed in the dark cloud TMC-1. In 
general we find good agreement between these observations and the abundances calculated 
in the molecular cloud model. 
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2.8. The disk model 



We have used the disk model of Id'Alessio et al.l (120011 ) with a mass accretion rate, M = 
10~ 8 M yr -1 to provide the physical parameters of the disk used here. The central star has 
a temperature of 4000 K, a mass of 0.7 M and a radius of 2.5 R . The surface density is 
25 gcm~ 2 at 10 AU and varies as l/R for R > 10 AU. The disk mass out to 400 AU is 0.063 
M . The grain size distribution is given by n(a) oc a -3 ' 5 with 0.005 /iin < a < 0.25 /zm. 
The density and temperature distributions are shown in Figure [H The gas temperature is 
set equal to the grain temperature, and the density and temperature are assumed to remain 
constant over the timescale considered in the chemical model. 



3. Results 

In the discussion that follows we look at the results of our models. Since the majority 
of molecules have been observed in the gas phase of disks we concentrate on these, but we 
also include a discussion of the grain chemistry and abundances because in some regions of 
the disk the solid phase is the dominant component, and because of the importance of grain 
chemistry in determining the gas phase composition and molecular deuteration. We look 
at how the molecular abundances vary with height above the midplane, and at the radial 
column density distribution. Finally, we compare our results with the available observations 
and use this to determine whether we can say anything about the nature of the desorption 
processes acting in the disk. 



3.1. The vertical molecular distribution 



3.1.1. Gas phase 



We begin by looking at the vertical abundance distributions at R = 250 AU and a time 
of 1 Myrs. In this section we concentrate on Model B; differences in the chemistry that 
arise because of desorption processes included in the other models are discussed in the next 
sections. 

The fractional abundances as a function of height, z, above the midplane are displayed 
in Figure El The three layer structure (see Figure El) fou nd by previous authors (lAikawa et al. 
20021 ; IWillacy fc Langerll2000l ; lAikawa fc Herbstlll999bl ) is clearly seen, with most molecules 
having low abundances in the midplane (due to freezeout at z < 50 AU) and in the surface 
layers (due to photodissociation at z > 100 AU), and abundance peaks in the molecular 
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layer, which for this model is between 50 and 95 AU above the midplane. 

The midplane is not completely devoid of molecules since desorption by CRH can main- 
tain low levels of CO and N 2 in this region. The reaction of these two species with Hjj~ and its 
deuterated isotopomers results in the formation of N2H + , N2D + , HCO + and DCO + . Since 
the temperature is cold, there is a high level of deuteration of producing high D/H ratios 
in the daughter molecules. The presence of CO and N 2 in the gas reduces the deuteration of 
H^" compared to Model A (which does not have non-thermal desorption) from Djj" /Hg ~ 70 
to a value of ~ 25. Deuterated Hjj~ forms from the reaction of HD with the previous molecule 
in the chain e.g. H 2 D + reacts with HD to form HD^ . If CO and N 2 are present in the gas, Hjj" 
and its deuterated isotopomers are more likely to react with t hese molecules than wi t h HD , 



reducing the formation rate of the deuterated molecules (see ICeccarelli &: Dominikl (120051 ) 
for details). 

The molecular layer begins at the point where thermal desorption can remove some 
molecules, mainly CO, N 2 and CH4 from the grains. The upper boundary is set by pho- 
todissociation. In the molecular layer most of the CO cycles between its gaseous and solid 
phases, but a small proportion is broken up by reaction with He + releasing C + and O, some 
of which go on to form molecules in the gas such as H 2 and H 2 CO. At z ~ 90 AU photodis- 
sociation can begin to destroy CO, forming carbon and oxygen atoms which do not reform 
CO but instead are converted into other molecules such as H 2 and hydrocarbons. These 
subsequently freezeout and are not sufficiently volatile to be thermally desorbed. Carbon 
and oxygen are thus removed from the gas, leading to a sharp decrease in the abundance of 
CO at z ~ 90 AU. At the very top of the disk most of the carbon is in its ionic form, since 
it is efficiently produced by photoionization and the temperature is warm enough that its 
freezeout is not efficient. 

N 2 exists in a narrower region than CO. CRH ensures that it is present in low abundances 
in the midplane. Above this, just below the level at which thermal desorption becomes 
efficient (~ 50 AU), its fractional abundance shows a decrease from its midplane value. The 
loss of N 2 is a result of its conversion into N 2 H + and N 2 D + . These ions can be destroyed either 
by reaction with CO or by recombination with electrons. At z ~ 45 AU, the main process is 
recombination, which produces mostly nitrogen atoms with NH and ND. These freezeout and 
are quickly hydrogenated to form ices of ammonia and its deuterated equivalents. Once CO 
begins to thermally desorb, its reaction with N 2 H + and N 2 D + becomes their main destruction 
route, producing N 2 , rather than N and NH or ND, and so less nitrogen is lost by conversion 
into non-volatile ices. 



Moving to higher z, there is a relatively narrow peak in the abundance of N 2 due to 
thermal desorption. This is followed by a rapid decrease in abundance at z ~ 65 AU where 
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photodissociation becomes efficient and the molecule is destroyed. The resulting nitrogen 
atoms can freezeout onto grains where they react quickly with hydrogen atoms to form NH, 
the first step towards forming ammonia ice. However NH is very volatile, with a binding 



energy of only 604 K (lAUen &: Robinson! 119771 ) and so can be easily returned to the gas (if 



the temperature is high enough) before hydrogenation can occur. The volatility of NH drives 
the nitrogen chemistry in the gas at z > 100 AU. Formation of gaseous molecules such as 
HCN and CN can proceed via 

NH CN+ S HCN+ ^ HCNH+ HCN CN (10) 

The photodissociation of CN produces nitrogen atoms and the cycle starts again. This 
process is sufficiently rapid to maintain low abundances of HCN in the gas phase, even in 
regions with high UV fields. 

Several other molecules also have small abundance peaks at the surface of the disk e.g. 
CH^ and H 2 (Figure [3]). These are formed in the warm gas by efficient reactions involving 
the atoms and ions produced by the photodissociation of molecules. The number of molecules 
produced is relatively small and does not greatly affect the calculated column densities. 

The effects of the increase in temperature with z on the D/H ratios can be seen in 
Figure [3j H3 is highly deuterated at z < 50 AU, with the most abundant isotopomer being 
Djj\ At higher z, Hjj~ is more abundant and of the deuterated forms only H 2 D + is important. 
The peak in Hjj" between z = 90 and 130 AU is due to the heavy depletion of molecules 
that would normally destroy it e.g. CO and N 2 . The variation in D/H ratio with z for Hjj~ 
translates into variations in the D/H ratio of other molecules e.g. DCO + /HCO + > 1 at z < 
50 AU, but at z > 50 AU HCO + is the dominant isotopomer. 



3.1.2. How photodesorption affects the calculated gas phase abundances and the molecular 

deuteration 

Figure H] shows the fractional abundances at 250 AU as calculated in Model C which 
includes both photodesorption and CRH. (Model D, which includes only photodesorption 
has a similar abundance distribution to Model C for z > 60 AU, but with much lower 
gas phase abundances in the midplane.) Photodesorption provides an efficient means of 
returning accreted species to the gas, especially important for those molecules which are 
too strongly bound to be affected by thermal desorption. It produces higher fractional 
abundances than Model B in the molecular layer, where the UV field is high enough that 
molecules can be removed from the grains but low enough that they can still survive in the 
gas. Photodesorption also increases the depth of the molecular layer. 
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At R = 250 AU, photodesorption affects the vertical distribution of most molecules, even 
those that can thermally desorb. Without photodesorption, elements are gradually removed 
from the gas by being incorporated into less volatile molecules e.g. carbon is removed from 
CO and converted into hydrocarbon ices, and oxygen and nitrogen form H 2 and NH 3 
ices respectively. When photodesorption is included we see large increases in the extent of 
the CO and N 2 layers. The most dramatic differences are for the less volatile species. For 
example, the peak fractional abundance of HCN increases to 2.5 x 10 -8 in Model C from 6 x 
10~ n in Model B. NH 3 and H 2 show dramatic increases in both peak fractional abundance 
and column density. All of these molecules can form rapidly on the grain surface, where 
deuteration is also efficient and so high column densities of HDO, D 2 0, NH 2 D etc. are also 
seen. Similarly, N(H 2 CO) is increased by grain formation and desorption, although to a 
lesser extent than some of the other molecules. 



3.1.3. The effectiveness of cosmic ray heating 



Desorption due to cosmic ray heating (CRH) was first suggeste d as a means of pr eventing 
the complete depletion of CO from the gas in cold, prestellar c ores (ILeger et al.lll985T). It was 
extended to cover other molecules i n molecular cloud models (IHasegawa & Herbstlll9 93l and 



has been included in disk models by lWillacy fc Langerl (120001 ) and I Willacy et al 



(l2006h . Here 



we examine whether its inclusion affects the calculated abundances and column densities. 

A comparison of the column densities calculated at R = 250 AU for Models A and B 
(Table H|) shows that there is little difference for non-deuterated molecules. However, several 
deuterated molecules are affected by CRH, because this allows molecules to be exist in the 
gas in the cold midplane of the disk where deuteration is efficient. For example, in Model 
B we see high abundances of DCO + at z < 50 AU (Figure EJ. In Model A, CO is not 
present in the gas in the midplane, and therefore neither is DCO + . At z > 50 AU, where the 
temperature is high enough for CO to thermally desorb, the temperature is also high enough 
that deuteration is less efficient and so little DCO + is present in Model A. A similar effect is 
seen for deuterated nitrogen-bearing molecules such as N 2 D + . This difference in D/H ratios 
provides a possible means of using the observations to differentiate between the models as 
shown below. 
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3.1.4- Vertical abundance distribution of the ices 

Figure [5] shows how the mantle abundances vary with z at R — 250 AU for Models B 
and C. In both models, the most volatile species have sharp decreases in abundance at z ~ 
60 - 70 AU where thermal desorption is able to return them to the gas. This corresponds to 
the peak in gas phase abundance seen in Figure [3J In Model B, H 2 CO and its isotopomers, 
which have a higher binding energy than CO, N 2 and CH4, are removed from the grains 
at z = 128 AU. In Model C, all species leave the grain between z = 80 - 100 AU due to 
photodesorption. 

In Model B many of the less volatile ice species show little change in abundance with 
height above the midplane, but there are some interesting variations, especially for molecules 
that contain two or more deuterium atoms. In the case of D2O, this molecule forms on grains 
by the addition of deuterium atoms to an accreted oxygen atom, but hydrogenation competes 
with the deuteration resulting in the production of H 2 and HDO as well as D 2 0. The 
proportion of oxygen atoms ending up in D 2 depends on the relative abundance of gaseous 
atomic hydrogen and deuterium. At z < D/H is high, ensuring that deuterium addition is 
efficient. At higher z, photodissociation begins to break down the main molecular reservoirs 
of hydrogen and deuterium, namely H 2 , HD and D 2 . Since there is more hydrogen tied up in 
these molecules than deuterium, their photodissociation results in a decrease in the atomic 
D/H ratio, and hence deuterium addition to accreted oxygen atoms becomes less competitive 
with hydrogen addition. The rate of formation of D 2 falls and we see this as a decrease 
in its abundance above z = 90 AU. The same processes affect the abundances of NHD 2 and 
ND 3 which are also formed on the grains. 

Several molecules e.g. H 2 0, HDO, HCN and DCN show peaks in abundance between 
z = 50 and 150 AU. In this region, CO is abundant in the gas phase, where it can be 
broken up into C + and oxygen atoms by reaction with He + . Some of the oxygen atoms 
freezeout and react with hydrogen or deuterium, producing abundance peaks in H 2 and 
HDO. Similarly, the freezeout of C + can result in its reaction with nitrogen atoms on the 
grain surface to form CN and subsequently HCN and DCN. There is also a peak in the 
abundance of hydrocarbons in the ice in this region. These form efficiently on the grains by 
the reaction of carbon (produced by the photodissociation of CO in the gas) and hydrogen. 
Their abundance falls at z > 140 AU, where the majority of carbon exists as gas phase ions. 

Model C shows similar behavior at z < 80 AU. The peak in C0 2 ice abundance just 
below this is due to its efficient formation in the gas by the reaction of OH and CO. The 
C0 2 then freezes out. Photodesorption is active at this z but not sufficiently so to remove 
all of the accreted mantles. Above z = 80 AU photodesorption dominates and the mantles 
are removed by z = 130 AU. 
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3.1.5. How assumptions about the binding energy of hydrogen atoms affect the results 
As discussed in Section 



600 K flCazaux fc Tielens 



3.1.51 we have elected to use a relatively high value of -Ed(H) 



20021 ). This is almost double our previously assumed value of 
350 K, and results in more efficient H 2 , HD and D 2 formation at temperatures above 15 K, 
with correspondingly lower gas phase abundances of atomic hydrogen and deuterium. The 
binding energy not only controls the residence time of hydrogen atoms on the grain, but 
also affects their ability to react. Hence the higher binding energy used here will increase 
the time that a hydrogen atom remains on the surface but will also reduce the rate at which 
hydrogenation processes can occur. 

To determine how Ed(H) affects the chemistry we reran Model B at R = 250 AU but this 
time with Ed(H) = 350 K and Ed(D) = 371 K. All other binding energies and assumptions 
remain the same. The reduction in Ed(H) increases the abundance of atomic hydrogen at 
z > 50 AU. For example, at z = 70 AU (within the molecular layer) s(H) increases from 
2.5 x 10~ 7 for E D (H) = 600 K to 4.7 x 10~ 4 for E D (H) = 350 K. A similar increase is seen 
in the atomic deuterium abundance. The column density of H doesn't change but that of D 
increases from 6.5 x 10 15 cm" 2 to 3.5 x 10 16 cm -2 . 

The changes in H and D abundance lead to changes in the abundances and column 
densities of other gaseous molecules. (The column densities of molecules in the ice mantles 
do not show much variation with the change in Ed(H)). In the gas phase, with Ed(H) =350 
K N(H 2 0) increases by - a factor of 10, as does N(H 2 CO), N(HDCO) and N(NH 2 ). The 
largest increase is for CH 3 OH whose column density increases from 6.5 x 10 5 cm -2 to 2.2 
x 10 7 cm -2 . The reduction in binding energy also changes some deuteration levels, most 
noticeably DCO + /HCO + increases to 3.4 (up from 0.9 in Model B). 

Using the lower value of £7d(H) in Model D changes the column densities even more. 
There are large reductions in column densities of factors between 5 and 10 for many of the 
deuterated molecules, including D 2 0, all of the deuterated isotopomers of CH 4 and NH 3 , 
DCN, and C 2 D. For example, N(D 2 0) is reduced from 7 x 10 12 cm -2 to 8.7 x 10 11 and 
N(DCN) from 1.2 x 10 12 cm" 2 to 1.2 x 10 11 cm" 2 . DCO+/HCO+ is also a factor of 4 higher 
than in Model D. 

The choice of H-atom binding energy can therefore have quite an effect on the calculated 
abundances. We have chosen the larger value since the results for a few species are in better 
agreement with the observations: in particular using Ed(H) = 350 K produces DCO + /HCO + 
ratios which are much higher than observed. 
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3.2. The radial distribution of molecules 

The radial distribution of molecular column densities, as calculated in Model B, are 
shown in Figure [61 N(CO) is fairly flat with radius for R > 100 AU, but shows a sharp 
rise at R = 50 AU, where thermal desorption of this molecule is efficient throughout the 
vertical extent of the disk. N 2 has a more gradual increase because of its slightly lower 
binding energy which means that it can desorb at larger radii. As discussed in the previous 
section, the presence of these volatile species in the gas can lead to the production of other 
molecules, which are not themselves thermally desorbed. For example, photodissociation of 
N 2 and CO can release the elements required to form molecules such as HCN and HCO + , 
and the column densities of these species increase towards the central star. 

N(Dg ) falls off with decreasing radius, because of the increase in temperature which 
inhibits its formation from Hj~. The relative levels of deuteration in some other molecules 
such as HCO + , H 2 0, HCN and NH 3 also decrease with decreasing R. CH 4 and H 2 CO do 
not show this decrease in deuteration, because they can form from CH 2 D + , which is less 
temperature sensitive due to the higher energy barrier for the conversion of CH 2 D + back 
into CHjj" (Ea = 370 K compared to Ea = 230 K for the conversion of H 2 D + into Hg ). Hence 
deuteration via reactions of CH 2 D + is more efficient in the warmer regions of the disk, than 
deuteration via H 2 D + . 

The column density of H 2 falls off as the radius decreases until R = 200 AU where it 
begins to rise again. The decrease at R = 50 AU is due to the high photon field. 

The radial distribution of column densities is also affected by photodesorption (Model C: 
Figure [7]), which increases the column densities of molecules at larger radii. Little difference 
is seen in the column densities of volatile species in the two models at R = 50 AU, but at 
larger radii photodesorption increases both N(CO) and N(N 2 ). In contrast, some molecules 
are little affected at large radii by the inclusion of photodesorption. For example, DCO + 
and HCO + have the same column densities in both models, as do N 2 H + , CH 4 and their 
isotopomers. 

The distributions of HCN, DCN and CN are all related, with CN being a photodisso- 
ciation product of the other two molecules. Their column densities in Model C are fairly 
constant with radius for R > 200 AU, but all three molecules show decreases in column 
density at smaller radii, caused by the increase of photodissociation rates closer to the star. 
The column densities in Model C are much higher than in Model B and the radial distribu- 
tions are different. In Model B, CN, HCN, and DCN all peak towards the center of the disk, 
whereas in Model C CN and HCN show a slight decrease. 
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3.3. The effect of disk chemistry on D/H ratios 

3.3.1. D/H ratios in the gas phase 



Aikawa fc Herbstl (jl999al ) found that chemical processing in the disk can affect the 
molecular D/H ratios and they therefore suggested that the ratios observed in comets do 
not just reflect the values set in the parent molecular cloud, but are instead a combination 
of processing in the cloud, during the infall phase and in the disk. Similarly, we find that 
gas phase D/H ratios are greatly changed by the disk chemistry. The situation for the 
grain mantles is a little more complicated, with disk chemistry having large effects on some 
molecules, and very little on others. This will be discussed below. 

Table [5] lists the calculated D/H ratios at R = 250 AU for each model and compares 
them to the input values as determined in the molecular cloud model. In the midplane very 
few molecules remain in the gas. Those that are present show a considerable increase in 
deuteration ratios compared to the input values. The effect is most marked for multiply 
deuterated molecules. Thus D^/Hg increases from 3.1 x 10~ 3 in the molecular cloud to 25 
in the midplane of Models B and C, due to the high molecular depletions in the disk. This 
compares to H 2 D + /H^ in the same models, which increases from an input value of 0.16 to 
~ 0.95. The increase in the midplane Dg /Hg ratio is even higher in Models A and D where 
the depletion is also higher due to the absence of CRH. 

The radial variation of the ratios of the column densities of deuterated to non-deuterated 



molec ules in Model B (Figured]) shows a similar distribution to those of lAikawa fc Herbst 



(120011 ). There is a decrease in D/H with decreasing radius for molecules such as NHg, H2O 
and HCO + whose deuteration depends on the isotopomers of Hjj~. The deuteration of H^ 
decreases towards the star because of the increase in temperature. Molecules such as H 2 CO 
and CH 4 whose deuteration depends on CH3 are less affected by the temperature increase, 
since the barrier to the reverse reaction of 

CH+ + HD ^ CH 2 D+ + H 2 (11) 



at 370 K is higher than that for H 2 D+ + H 2 ^ H+ + HD (220 K). The ratio CH 2 D+/CH+ 
decreases towards the center of the disk, because there is a fall in HD in the molecular layer, 
caused partly by the increase in photodissociation as the UV field increases and partly by 
the efficient incorporation of deuterium into water and ammonia ices. 

With the inclusion of photodesorption the picture changes. This mechanism can alter the 
molecular D/H ratios by injecting molecules formed on the grains into the gas and therefore 
the gaseous D/H ratios partly depend on the ratios in the ice mantles. Some molecules 
do not show much variation in the D/H ratio with photodesorption e.g. HDCO/H 2 CO and 
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DCO + /HCO + are roughly the same in both Models B and C. 

In Models C and D the high abundance of molecules in the molecular layer keep the 
abundances of H^" and Djj~ in this region low. Hence most of the Hjj~ and Djj" are in the 
midplane where the temperatures are low and the deuteration is high leading to high D^/H^ 
ratios for these models. In Models A and B the higher depletion at z > 50 AU means that the 
destruction rate of Hjj~ is lower and we see a secondary peak in abundance for this molecule 
at z ~ 70 - 120 AU (Figure EJ. Consequently N(H^~) is higher and N(D^)/N(H^) is lower 
in these models. 

In Mo del B we find that DCN /HCN decreases slightly with decreasing radius in agree- 



ment with lAikawa fc Herbstl (120011 ). although our ratios are slightly higher than theirs. At 
large radii DCN forms from DCNH + which is produced by the reaction of either DCO + or 
with HNC. Hence DCN/HCN is dependent on the deuteration of Hj. At R = 50 AU, neutral- 
neutral reactions become important with DCN forming from D + HNC — > DCN + H and 
therefore the DCN/HCN ratio is dependent on the atomic D/H ratio. Another important 
route to both DCN and HCN at small radii is the reaction of N with either CH 2 D or CH 3 . 
DCN/HCN is also relatively constant with R for Model C, where both molecules form on 
the grain surfaces and then photodesorb. In this case their ratio is therefore dependent on 
the gas phase atomic D/H. 

In summary, chemical processing in the disk can affect the D/H ratios of some molecules. 
The effects is strongest for multiply deuterated molecules, and for ammonia. The radial 
variation of molecular deuteration depends on the formation route of the deuterated molecule, 
with species that are deuterated via reactions with H2D + , HD^~ and Djj~ showing a decrease 
in the deuteration with decreasing radius, in contrast to molecules deuterated by CH 2 D + , 
CHD^ and CDg . If photodesorption is included then the formation of deuterated molecules 
on the grains becomes more important. 



3.3.2. D/H ratios in the ice mantles 

In general, chemical processing in the disk has less effect on the D/H ratios of the ices 
than on the gas phase molecules. Looking first at the midplane (Table [5]) we find that the 
ratios calculated for the ices in Model A are close to the input values, with the exception of 
DCN: gr/HCN: gr. Model A also shows no change in the fractional abundances of most ice 
molecules. At the start of the model, most molecules have much higher abundances on the 
grains than in the gas, and so the freezeout of the gaseous molecules does not greatly affect 
the grain mantle abundances. The exceptions are for CO, N 2 , HCN and DCN. For CO and 
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N2 the initial gas phase abundances are high. In the case of HCN and DCN these molecules 
form on the grains from the freezeout and hydrogenation/deuteration of CN. 

CRH causes some changes in the ice deuteration in the midplane in Models B and 
C. As in the gas we find that the biggest changes are for multiply deuterated molecules. 
The break up of desorbed CO and N2 provides carbon ions and oxygen and nitrogen atoms 
which can freezeout and react on the grain surface, most likely with hydrogen or deuterium. 
The relative efficiency of these two processes depends on the atomic D/H ratio, which in 
the midplane is very high. Hence multiply deuterated ices are easily formed in this region. 
Ammonia ice is particularly highly deuterated because a lot of it is produced in the disk 
(due to the high input abundances of both gaseous and solid N 2 which provide a plentiful 
source of nitrogen atoms), whereas for H 2 and CH 4 the majority of these ices are produced 
in the molecular cloud phase and the disk chemistry makes a much lower contribution to 
both their abundances and their deuteration. 

If we calculate the molecular D/H ratios from the column densities of the ices, we 
again find that Models A and D show the least change from the molecular cloud values. In 
Models B and C the column density ratios show similar behavior to the midplane ratios, 
with multiply deuterated species and the isotopomers of ammonia being most affected. 

In summary we find that the inclusion of a non-thermal desorption process, such as 
CRH, which is effective in the midplane will affect the D/H ratios of some of the ice molecules 
and that this effect increases for multiply deuterated molecules. Such a process would be 
expected to alter the D/H ratios of the material that is later incorporated into Solar System 
bodies such as comets. In the absence of such a non-thermal desorption process, the ice 
mantles retain the D/H ratios set in the molecular cloud and are unaffected by the disk 
chemistry. The exception to this is ammonia, which is predicted to be highly deuterated in 
the ice by all of our models. At smaller radii, where thermal desorption is efficient, we would 
expect that disk chemistry will play a more important role in determining the D/H ratios of 
the ices. 

There are limitations to our model. We have not incorporated chemical processing of 
the grain mantles by cosmic rays or UV photons. UV will not be important in the heavily 
shielded midplane and so can be safely ignored in this region, but cosmic rays can penetrate 
and may have a chemical effect. Whether or not this can affect the D/H ratios is unclear. 
A further complication ignored here is processing of the grain mantles during the collapse 
phase as the molecular cloud material is accreted into the disk. An accretion shock forms at 



Lunine et al. 


( 


1991 





water ice coatings if they are infalling at radii of less than 30 AU, with the effect decreasing 
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with increasing radius (ILunine et al.lll99ll ). The grains will regain their ice mantles once they 
have passed through the shock, but processing in the shock could alter the composition and 
deuteration levels of the ices. Both these papers suggest that at least some of the material in 
comets comes directly from the inte r stella r medium. At the r adii u nder consideration here, it 
seems likely based on lLunine et al.l (1199 lh : IChick fc Cassenl (119971 ) that at least some of the 
grains will retain their volatiles as they are accreted into the disk. At present it is unclear 
what effect the accretion shock would have on the D/H ratios of the ice mantles on grains 
at smaller radii, where these mantles may have been desorbed in the shock and the released 
molecules processed in the gas before re-accreting onto the grains. 

It has been suggested that comets are composed of interstellar grains that have u nder- 
gone relatively little processing in the Solar Nebula (see review by llrvine et al.l |2000| . and 
references therein). One reason for supposing this is the similarity between the chemistry of 
comets and interstellar clouds, and in particular the similarity between the molecular D/H 
ratios. Dynamical studi es of comets have s hown that short period comets formed in the 
trans-Neptunian region (IDuncan et al.lll988l ) whereas long period comets formed closer to 
the Sun and then had their orbits perturbed by interaction with the giant planets, to eject 
them outwards to the Oort cloud. Our model is of the outer disk, with an inner radius well 
outside the region where long period comets - the only ones for which we have observed D/H 
ratios - are likely to have formed. Still we include the comet data here for comparison, while 
recognizing that models of the inner disk are requi red to properly cove r the comet formation 
region. In Hale-Bo pp, DCN/HCN = 2.3 x 10" 3 JMeier et alJll998ah. HDO/H^O = 3.3 x 
10~ 4 in Hale-Bopp (jMeier et al.lll998al ) and 2.9 x 10~ 4 in Hyakutake (IBockelee-Morvan et al. 



19981 ). Our model results are close to the input values with DCN: gr/HCN: gr ~ 6.6 x 10 3 
and HDO: gr /H 2 0: gr ~ 2.7 x 1(T 2 in the midplane at R = 50 AU (Model D). The 
DCN: gr/HCN: gr ratio is consistent with the observations but HDO: gr /H 2 0: gr is much 
higher than the observed value. The discrepancy could be a result of too efficient deuteration 
on grain surfaces, or it could be because our model considers a region further out in the disk 
than the comet formation radius. Processing at a small radii (with higher temperatures) 
could reduce the calculated HDO/H 2 value. 



4. Comparison to observations 

One reason for modeling the outer parts of the protoplanetary disk is that most current 
observations in the millimeter and submillimeter probe this region and therefore can provide 
a test of the models. Three disks in which several molecules (including deuterated forms) 
have been observed are LkCal5, DM Tau and TW Hya. 
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TW Hya: This object is the closest classical T Tauri star known. Its disk is almost face 



It has a mass accretion rate, M — 10 9 - 10 8 M yr 
disk mass is estimated from continuum observations to be 3 x 10 



2000). 



Kastner et al. 


20021. The 


3 x 10" 2 M r 


Wilner et al. 



• LkCalS: This s ource is one of t h e stro ngest T Tauri emitters in the millimeter found 
in the survey of iBeckwith et al. I dl990h. It is found located in Taurus and has an age 



Thi et a 



2(iQ^. The disk mass is - 0.03 M . 



is 650 AU (ISimon et al. 



temperature is 4365 K (IMuzerolle et al.l 120001 ) 



200 01) . The stellar m ass is 1 M ( ISimon et al.ll2000l ) and its 



of ~ 3-10 Myrs (ISimon et all 120001 ; 

and the mass accretion rate is 10~ 9 M Q yr -1 (IHartmann et al.lll998f). The dis k radius 



DM Tau: This object is also found in Ta urus. It has an age of ~ 5 Myrs. Its 
stellar mass is 0.55 (pinion et al.l 120001 ) and its stellar te mperature is 3720 K 
(IGuilloteau fc Dutreylll998l ). The disk has a m ass of 0.03 M (IGuilloteau fc Dutrey 
19981 ) , a radius of ~ 800 AU (ISimon et al.l 120001 ) and a mass accretion rate of 10~ 8 M 



yr 1 ( IHartmann et al.lll998l ). 



We compare our results with the available observations of both deuterated and non- 
deuterated species. The column densities derived from observations are very sensitive to 
assumptions made about the source. Single disk observations cannot resolve the disks and 
assumptions about the density and temperature structure of the dis ks must be made in 
order to determine abundances and column densities. For DM Tau, iDutrey et al.l (119971 ) 
determined the gas density distribution of a geometrically thin disk in hydrostatic equilibrium 
and used this to derive the average fractional abundances with res pect to H ? , as s uming that 
the fractional abundances were the same every where in the disk. lAikawa et al.l (120021 ) used 
the model and data from IDutrey et al.l (119971 ) and integrated vertically to determine the 
column densit ies qu oted in Table [61 For LkCal5 both single dish and interferometric data is 
available. |QJ (120011 ) used the Owens Valley Radio Interferometer (OVRO) array to observe 
this source and derived beam average column densities. The resolution of the array f~ 300 



AU at the distance of LkCal5) is such that the source is just resolved. iThi et al.l (120041 ) 
present single dish observations of LkCal5. They derive the column densities assuming that 
the radius of the disk is 450 AU. The disk radius can have a large effect on the calculated 
column densities with a change from Ri to R2 scaling the column densities by (i?i/i?2) 2 - 
The single dish column densities are significantly higher than those from the interferometer 
data due to the differences in the assumptions used to determine these numbers. 



Several deuterated molecules are now claimed to have been detected in protoplane- 
tary disks. However, it should be noted that some of these detections have recently been 
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(jCeccareUi et al 



disputed. iGuillqteau et al.l (120061 ) re-analyzed the spectra of HDO and H 2 D + in DM Tau 



20051 ) and found only a 2-a detection of H 2 D + and no evidence of HDO 



in this disk. They attributed the discrepancy between their work and that of the original 
authors to the use of an over-estimation of the continuum flux and the wrong systematic 
velocity by the latter. In the following discussion we have elected to use all the available 
data on deuterium molecules in disks for comparison with our models, including that which 
is still under dispute. 

The observations referred to here are of the outer disk, where the radius is greater than 
50 - 100 AU. Here we use our results at R = 250 AU as representative of the outer disk 
model and compare these to the observations. The results for Models B, C and D, together 
with the observations are given in Table O Good agreement is defined as less than a factor of 
5 difference between the models and the observations. We have not included Model A in this 
table since the column densities for this model are consistently lower than the others and 
much lower than observed, indicating that non-thermal desorption is required to account for 
the observations. 

For each source we find that Models C and D give better agreement with the observations 
than Model B. In particular, strongly bound molecules such as HDO require the presence of 
photodesorption in order to remain in the gas phase. 



Ceccarelli et al.l (120051 ) observed HDO in absorption in DM Tau using the JCMT and 
determined N(HDO) = 1.6 x 10 13 cm" 2 . In Model B, N(HDO) is too low at 9.3 x 10 9 
cm -2 to account for the observations, but in Models C and D photodesorption can keep a 
much higher abundance of HDO in the gas. These models have N(HDO) = 2.0 x 10 14 cm -2 , 
somewhat higher than the observed value, but demonstrating that photodesorption is an 
efficient means of returning HDO (and other molecules) to the gas. 



Dominik et al.l (120051 ) also looked the effects of photodesorption on water and its iso- 



topomers. Using a simplified model they calculated the column density of H 2 averaged 
across a face on disk to be 1.6 x 10 15 cm -2 , in exce llent agreement with ou r model. Com- 
bining their model results with the observations of ICeccarelli et al.l (120051 ). they estimate 
HDO/H 2 ~ 0.01. Our models calculate much higher ratios with HDO/H 2 = 0.13 in 
Models C and D. These calculated ratios are also much higher than are observed in Solar 
System objects and may point to deuteration on the grains being less efficient than assumed 
here. 



There is another reported detection of HDO in a protostellar disk. |OJ (120011 ) observed 



HDO in LkCal5 using the OVRO interferometer and found N(HDO) = 2-7 xlO 14 cm" 2 , 
with an intensity peak that was offset from the central star. As in DM Tau, models that 
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include photodesorption give the best agreement with the observations. 



Qj (120011 ) also observed DCN in LkCal5 with N(DCN) 
0.01. Model B calculates a low value of N(DCN) = 8.9 xlO 8 cm" 2 



10 13 cm" 2 and DCN/HCN 



but Models C and D have 



N(DCN) = 1.2 xlO 12 cm" 2 (DCN/HCN = 0.036), in good agreement with the observations. 



H 2 D + has been observed in DM Tau (ICeccarelli et al.ll2004l ) with a column density of 
8.8 xlO 12 cm -2 , and a midplane fractional abundance x(H 2 D + ) = 3.4 xl0~ 10 . In all our 
models we find x(H 2 D+) ~ few xlO" 12 in the midplane for R > 100 AU, and N(H 2 D+) ~ 
few xlO 11 cm -2 . Our calculated N(H 2 D + ) is therefore a factor of 10 lower than observed 
and our midplane fractional abundance more than a factor of 100 lower. 

DC O+ and HCO+ haye bot h been detected in two sources: DM Ta u (DCO+/HCO+ 



4 xlQ- 3 : iGuilloteau et al]l2006h and TW Hya (DCO+/HCO+ = 0.035: Ivan Dishoeck et al. 



20031 ). The formation of these molecular ions depend s on the presence of CO in the gas. The 



column density of C O in DM Tau is 5. 7 xlO 1 6 cm 2 (IDutrey et al.lll997l ) in good agreement 



with all our models. IGuilloteau et al.l (120061 ) do not give a column density for DCO + , but 
in TW Hya N(DCO + ) = 3 xlO 11 cm -2 compared to the calculated values of ~ 3.1 xlO 12 in 
Models B and C and 4.8 xlO 11 in Model D. In the same source, N(HCO + ) = 8.5 xlO 12 cm" 2 
and the model values, which range from 2.8 - 3.3 xlO 12 cm -2 , are in good agreement with 
this. In our models the desorption process that affects the calculated ratios most is CRH. 
This allows CO to be present in the cold midplane where deuteration is most efficient. When 
CRH is included N(DCO+) is high and our calculated DCO+/HCO+ ~ 1 (Models B and C). 
Without CRH (Model D) both N(DCO + ) and the ratio of 0.18 are in agreement with the 
TW Hya observations. However for this model DCO + /HCO + is 42 times higher than seen 
in DM Tau. The observations and modeling of the DCO + /HCO + ratio suggests that CO is 
not present in the midplane, but instead does not begin to become abundant until the disk 
is warm enough that deuteration is not so efficient, as happens in Model D. Based on this 
we suggest that CRH is not acting in disks, or if it is, it is less efficient than assumed in our 
(and other) models. 



Qi et all (|2003h observed N(N 2 H H 



3 cm 2 in LkCal5. The highest column 



= 3.1 x 10 

density found in any of our models is 2.8 x 10 10 (Model C), ~ 1100 times lower than 
observed. The discrepancy is even larger for Model D, where N(N 2 H + ) = 1.1 xlO 10 cm~ 2 . 
The increase in N(N 2 H + ) in models that include CRH suggests that a high column density 
of N 2 H + requires that N 2 is present in the midplane of the disk. 

So from a comparison of observations with our models we can see that we require N 2 to 
be present in the midplane, but CO to be absent. In the current models the binding energies 
of these two molecules are similar and therefore they have similar desorption rates. In order 
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to account for the observations we would need j?p(N 2 ) to be somewhat smaller than EpiCO). 
The values we use are taken from lOberg et al.l (120051 ) who measured the binding energy of 
N 2 on N 2 ice, CO on CO ice and the binding energies of both molecules on a mixture of CO 
and N 2 ices. However, in the interstellar medium and protostellar disks the main component 
of the ice mantles is water and it is possible th at the binding energ ies of N 2 and CO to water 
ice are very different from those measured by lOberg et al.l (120051 ) . Various theoretical and 
experimental work points towards this being the case. For example, iHerbst fc Cuppenl (120061 ) 
take En (CO) and £ , p(N 2 ) to be 11 50 K and 100 K re spectively, based o n the laboratory 
work of ICollings et al.l (120031 ) and lAyotte et al.l (120011 ). Calculations by lAl-Halabi et al. 



(120041 ) find -Ed (CO) to be 1010 K on amorph ous ice and 1215 K on crystalline ice. These 
values are rather more than the lOberg et al.l values. Increasing the binding energies will 
reduce the efficiency of both thermal desorption and CRH (for which rates are calculated 
using E D ). Also increasing the relative value of E D {CO) compared to E D (N 2 ) will mean 
that N 2 could be desorbed in regions where CO is retained on the grains. For example, using 
E D (CO) = 1200 K and E D (N 2 ) = 1000 K, we find midplane abundances x(N 2 ) = 3.2 xl0~ 10 
and x(CO ) = 5.2 xlO" 10 at R = 250 AU, compared to 1.4 xl0~ 9 and 6.0 xlO" 8 respectively 



using the lOberg et al.l binding energies. In terms of column densities, the higher binding 
energies have N(N 2 H + ) = 6.2 xlO 10 cm" 2 (an increase of a factor of 2.4 over the value in 
Model C), while N(CO) = 1.4 xlO 17 cm" 2 and DCO+/HCO+ = 0.23. Thus with the higher 
binding energies, N(CO) is still in agreement with the observations and N(N 2 H + ) increases 
but is still a factor of ~ 500 lower than the observed value in LkCal5. A higher column 
density could be achieved if N 2 is more efficiently desorbed in the midplane. 

Even with the increased CO binding energy, DCO + /HCO + is still much higher than 
observed, because CRH can still desorb some CO in the midplane. Another possibility is 
that CO does not desorb by CRH a s efficiently as we assu me here. Our rates a re bas ed on the 
work of IHasegawa fc Herbstl (Il993l ). but recent work by iBringa fc Johnson! (120041 ) suggests 
that theses rates maybe an over-estimate for CO and an underestimate for H 2 0. According 
to Bringa and Johnson the d esorption rate for CO is a factor of 10 lower than estimated 
by IHasegawa fc Herbstl (119931 ) which would reduce the abundance of CO in the midplane 
considerably. There is no data for N 2 , but if it is able to desorb at the Hasegawa & Herbst 
rate, then this would allow there to be a midplane layer of N 2 that is also lacking in CO. 

On the basis of our modeling we find that efficient non-thermal desorption such as pho- 
todesorption must be active in the surface and molecular layers to account for the observed 
column densities of molecules such as HDO and HCN which would otherwise be expected 
to be accreted onto the grains. The observed DCO + /HCO + ratio indicates that CO is not 
present in the midplane where deuteration is efficient, ruling out desorption of CO in this 
region that is as efficient as the CRH included here. However, some desorption must be 
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occurring in the midplane to account for observations of N2H + . Overall Model D gives the 
best agreement with observations, but in order to account for N 2 H + , N 2 must be able to 
desorb in the midplane. 



5. Summary 

We have presented results from new models of the deuterium chemistry in disks which 
include grain mantle chemistry and multiply deuterated molecules, and have considered how 
desorption processes affect the calculated abundances and the molecular deuteration. Our 
main findings are as follow: 

1. Very high levels of deuteration can be achieved in the midplane, where freezeout re- 
moves most molecules from the gas. 

2. Efficient desorption is required in the molecular layer to acco unt for the observat i ons o f 



mol ecules such as HDO. I n common with other authors e.g. IWillacy &: Langerl (120001 ) 



and iDominik et al.l (120051 ). we find that photodesorption is able to account for the 



observed column densities. 

3. Chemical processing in the disk, including grain surface chemistry, can make great 
changes to the deuteration of gas phase molecules compared to the input molecular 
cloud values. This is partly a consequence of processing in the gas and partly due 
to the formation of molecules on the grains and their subsequent desorption. We find 
that DCO + /HCO + is very sensitive to the desorption processes included in the models. 
From a comparison of observational data for this molecule and the models we find that 
CO cannot be present in the midplane, suggesting either that the rate of cosmic ray 
desorption of this molecule is less than assumed in these models, or that its binding 
energy on water ice is much higher than the value adopted here. 

4. From comparison of the model results in the midplane at 50 AU we find that some 
molecules found in comets have D/H ratios which reflect the input molecular cloud 
abundances, but that others, notably NH 3 have undergone considerable chemical pro- 
cessing in the disk. Based on this it seems that comets could be comprised of some 
molecules that have their origin in the parent molecular cloud and others that were 
formed in the disk. 



The deuterium chemistry provides a means of testing the models against observations 
and of distinguishing among the various processes that might be acting in the disk. As 
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observations continue to be acquired and gain in sensitivity and resolution, the ability to use 
models and observations together to understand the chemistry of protoplanetary disks will 
be greatly enhanced. 

This research was conducted at the Jet Propulsion Laboratory, California Institute of 
Technology under contract with the National Aeronautics and Space Administration. Partial 
support was provided by the NASA TPF Foundation Science Program. 
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Table 1. The binding energies {Ed) used to determine the thermal desorption rates of the 

abundant mantle species. 



Species 


E D (K) 


Reference 


Species 


E D (K) 


Reference 


H 


600 


1 


D 


621 


2 


H 2 


315 


8 


C 


800 


3 


CH 


645 
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CH 2 


956 


9 


CH 3 


1158 


9 


CO 


855 


4 


co 2 


2860 


5 


H 2 CO 


1760 
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CH3OH 


4240 


7 





800 


3 


o 2 
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3 


OH 


1259 


9 


H 2 


4820 


6 


N 


800 


3 


N 2 


790 


4 


NH 
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9 


NH 2 
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9 


NH 3 


3080 
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Table 2. Elemental abundances used in the model. Abundances are given with respect to 
the total abundance of hydrogen atoms, nu = 2 n(H 2 ) + n(H). Initially all deuterium is 

contained in HD, and all hydrogen in H 2 . 



Element 


Abundance 




n(x)/riH 


H 2 


0.5 


He 


0.14 


HD 


1.6 x 10~ 5 





1.76 x 10~ 4 


C+ 


7.3 x 10~ 5 


N 


2.14 xlO" 5 


Fe 


3.0 xlO" 9 


Mg 


7.0 xlO" 9 
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Table 3. Input abundances for the disk model as determined by a molecular cloud model 
at 1 Myrs. The abundances are given relative to the total abundance of hydrogen % = 2 
n(B.2) + ra(H). For comparison we have included the abundances observed at the 
cyanopolyyne peak in the molecular cloud TMC-1. 



Molecule Fractional abundance TMC-1 
Gas Grain Gas phase 
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3.6 


(-9) 


2.4 (-6) 






D 2 


1.6 


(-11) 


1.8 (-8) 






N 


4.5 


(-7) 








N 2 


2.2 


(-6) 


7.3 (-8) 






NO 


1.1 


(-7) 




3.0 


(-8) 1 


CN 


2.5 


(-8) 




5.0 


(-9) 4 


HCN 


1.5 


(-8) 


4.6 (-6) 


2.0 


(-8) 1 


DCN 


4.9 


(-10) 


2.8 (-8) 


2.2 


(-10) 2 < 3 


HNC 


1.0 


(-8) 


7.2 (-7) 


2.0 


(-8) 1 


DNC 


1.2 


(-10) 


3.0 (-9) 


3.0 


(.10)3,5 


C 3 N 


2.9 


(-9) 




6.0 


(-10) 9 

(-8) 4 


HC 3 N 


5.8 


(-8) 


3.7 (-7) 


2.0 



Table 3 — Continued 



1V1U1LL U1C 


Fractional abundance 




TMC-1 




Gas 


Grain 




VJCIO jJIlCliOL. 


DC 3 N 


1.3 (-9) 


8.0 (-9) 


6.0 


(-10) - 2.0 (-9) 6 


CHgCN 


1.0 (-9) 


1.5 (-8) 


6.0 


(-10) 9 


NH 3 




l.U (-t>) 


2.0 


(-8)^ 4 


NH 2 D 


4.4 (-11) 


2.4 (-7) 


2.2 


(-10) 2 


NHD 2 


5.6 (-13) 


9.2 (-10) 






ND 3 


9.2 (-14) 


1.4 (-11) 






N 2 H+ 


2.8 (-10) 




4.0 


(-10) 1 


N 2 D+ 


2.2 (-11) 




3.2 


(-H) 2 


C 


7.8 (-7) 








c+ 


1.8 (-8) 








CH 


5.2 (-9) 




2.0 


(-8) 1 


C 2 H 


1.2 (-9) 




2.0 


(-8) 7 


C 2 D 


3.5 (-11) 




2.0 


(-10) 8 


C 3 H 


5.9 (-7) 




1.0 


(-8) 7 


C 2 H 2 


1.7 (-8) 


1.9 (-6) 






CH 4 


2.7 (-7) 


1.4 (-6) 






CH 3 D 


1.7 (-8) 


9.2 (-8) 






CH 2 D 2 


9.3 (-10) 


2.6 (-9) 






CHD 3 


3.3 (-11) 


1.8 (-10) 






CD 4 


1.2 (-13) 


6.4 (-13) 
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Table 4. The column densities, N(x), calculated at R — 250 AU and a time of 1 Myrs 
(a(b) represents a x 10 6 cm' 2 ). All models include thermal desorption and this is the only 
desorption process included in Model A. Model B also includes CRH, Model C includes 
CRH and photodesorption and Model D only includes photodesorption. The column 
density of H 2 = 4.8 x 10 22 cm -2 for all models. 



Molecule N(x) 

Model A Model B Model C Model D 



Gas species 



HD 


5.2 


(17) 


2.9 


(17) 


2.8 


(17) 


5.1 


(17) 


CO 


1.2 


(17) 


1.2 


(17) 


2.9 


(17) 


2.7 


(17) 


H 2 


2.7 


(10) 


4.0 


(10) 


1.6 


(15) 


1.6 


(15) 

yj-'-j ) 


HDO 


1.4 


(9) 


9.3 


(9) 


2.0 


(14) 


2.0 


(14) 


D 2 


1.5 


(7) 


1.3 


(8) 


7.1 


(12) 


7.0 


(12) 


co 2 


9.8 


(9) 


1.0 


(10) 


1.1 


(14) 


1.1 


(14) 


CH 4 


2.9 


(14) 


3.3 


(14) 


2.9 


(14) 


2.5 


(14) 


CH 3 D 


3.7 


(13) 


4.5 


(13) 


3.4 


(13) 


2.7 


(13) 


CH 2 D 2 


3.0 


(12) 


3.8 


(12) 


2.5 


(12) 


1.7 


(12) 


CHD 3 


1.1 


(11) 


1.8 


(11) 

V ) 


1.1 


(11) 

V ) 


4.6 


(10) 


CD 4 


1.3 


(9) 


6.9 


(9) 


5.9 


(9) 


2.6 


(8) 


HNC 


2.0 


(10) 


3.7 


(10) 


3.2 


(12) 


3.1 


(12) 


DNC 


8.9 


(7) 


2.2 


(8) 


8.5 


(10) 


8.6 


(10) 


HCN 


1.7 


(10) 


1.8 


(10) 


3.3 


(13) 


3.3 


(13) 


DCN 


2.3 


(8) 


8.9 


(8) 


1.2 


(12) 


1.2 


(12) 


CN 


3.7 


(11) 


3.8 


(11) 

V ) 


3.6 


(13) 


3.6 


(13) 


H 2 CO 


3.3 


(11) 


3.3 


(11) 


3.9 


(12) 


3.9 


(12) 


HDCO 


6.7 


(9) 


7.5 


(9) 


8.6 


(10) 


8.5 


(10) 


D 2 CO 


6.0 


(8) 


6.5 


(8) 


1.2 


(9) 


1.2 


(9) 


HC 3 N 


1.5 


(9) 


1.6 


(9) 


7.0 


(11) 


7.0 


(11) 


DC 3 N 


4.6 


(7) 


4.9 


(7) 


2.1 


(10) 


2.2 


(10) 


N 2 


3.8 


(15) 


3.9 


(15) 


2.2 


(16) 


2.2 


(16) 


N 2 H+ 


6.3 


(9) 


2.3 


(10) 


2.8 


(10) 


1.1 


(10) 


N 2 D+ 


3.2 


(9) 


1.6 


(11) 


1.6 


(11) 


3.8 


(9) 


CH 3 CN 


2.1 


(7) 


2.1 


(7) 


4.9 


(10) 


4.8 


(10) 


CH 3 OH 


6.4 


(5) 


6.5 


(5) 


1.5 


(11) 


1.5 


(11) 


CH 3 OD 


7.9 


(3) 


1.2 


(4) 


1.2 


(10) 


1.2 


(10) 


CH 2 DOH 


3.3 


(4) 


3.5 


(4) 


1.3 


(10) 


1.3 


(10) 


H 3 + 


2.8 


(12) 


3.1 


(12) 


8.9 


(11) 


5.3 


(11) 


H 2 D+ 


5.3 


(11) 


8.2 


(11) 


7.7 


(11) 


4.8 


(11) 


HD+ 


6.6 


(11) 


1.1 


(12) 


1.1 


(12) 


6.6 


(11) 


D 3 + 


2.6 


(13) 


1.3 


(13) 


1.3 


(13) 


2.6 


(13) 


HCO+ 


2.9 


(12) 


3.3 


(12) 


3.2 


(12) 


2.8 


(12) 


DCO+ 


3.4 


(11) 


3.0 


(12) 


3.2 


(12) 


4.8 


(11) 


CH+ 


9.2 


(10) 


9.2 


(10) 


1.1 


(11) 


1.1 


(11) 
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Table 4 — Continued 



Molecule 



Model A 



N(x) 

Model B Model C 



Model D 



V I 1 ! W 


1 n 


I 9 ] 


1 n 

l.U 


(y) 


A 9 
4.Z 


(?) 


a n 

4.U 


1°; 




A A 
4.4 


(') 


A Q 
4.0 


{') 


1 n 
1 .u 


I7'\ 
{') 


1 

y. 1 


(6) 


pn+ 


y.o 


(K\ 
( 5 ) 


y.o 


(5) 


1 7 
1 . 1 


(5) 


1.0 


(&J 


O2W 


0.0 


(12) 


Q ft 
O.O 


(12) 


7 9 
/ .z 


(12) 


7 9 
( .z 


(12) 


02-D 


o.U 


ci ni 
(1U) 


Q 1 
O.l 


(10) 


9 9 

z.z 




9 1 

z. 1 


I 11 ) 


INM3 


7 
Z. f 






z.y 




1 A 
1 .4 


CI 

(13) 


1 A 
1.4 


(13) 


IN X12-L' 


0.0 




1 Q 

i.y 




1 .0 


(1^2) 


1 9 
1 .Z 


[iZ) 


INXIU2 


z.o 




1 k 
1.0 




9 A 
Z.O 


(1U) 


9 Q 
Z.o 


(1U) 


INJJ3 


1.0 




A 
4.y 


(°) 


/I 9 
4.Z 


(9) 


Q 9 
O.Z 


(9) 


| I 

XI 


Q Q 
O.O 


k iy j 


4.0 


k iy ) 


/I Q 


fi ol 
(iy ) 


/I Q 
4.0 


(iy ) 


U 


D.O 


(lt>) 


0.0 


(lb) 


1 .0 


(10) 


1 3 
l.O 


(lb) 


IVlantle species 




















/I 1 
4.1 


(~\ ai 
(18) 


a n 
4.U 


/"1 Ql 

(18) 


4.U 


(18) 


/I 1 

4. 1 


(18) 


XI2U 


1 9 
1 .Z 


(iyj 


1 9 
1 .Z 


^iy j 


1 9 
1 .Z 


^iy j 


1 9 
1 .Z 


^iy j 


HDO 


2.7 


(17) 


3.2 


(17) 


3.3 


(17) 


2.8 


(17) 


D 2 


4.0 


(15) 


1.9 


(16) 


1.9 


(16) 


4.1 


(15) 


CH 4 


2.5 


(17) 


2.7 


(17) 


2.7 


(17) 


2.4 


(17) 


CH 3 D 


2.2 


(16) 


5.5 


(16) 


5.4 


(16) 


2.1 


(16) 


CH 2 D 2 


1.9 


(15) 


1.1 


(16) 


1.1 


(16) 


1.8 


(15) 


HCN 


4.5 


(17) 


4.9 


(17) 


4.7 


(17) 


4.7 


(17) 


DCN 


9.5 


(15) 


8.9 


(15) 


8.3 


(15) 


8.9 


(15) 


H 2 CO 


3.9 


(17) 


3.9 


(17) 


3.8 


(17) 


3.7 


(17) 


HDCO 


3.3 


(15) 


3.8 


(15) 


3.6 


(15) 


3.0 


(15) 


D 2 CO 


4.1 


(13) 


4.2 


(13) 


3.5 


(13) 


3.6 


(13) 


NH 3 


1.0 


(18) 


1.0 


(18) 


1.0 


(18) 


9.6 


(17) 


NH 2 D 


3.4 


(16) 


2.3 


(17) 


2.3 


(17) 


3.3 


(16) 


NHD 2 


2.1 


(15) 


5.2 


(16) 


5.2 


(16) 


2.0 


(15) 


CH3OH 


4.3 


(15) 


4.3 


(15) 


3.9 


(15) 


3.9 


(15) 


CH 2 DOH 


3.1 


(14) 


3.2 


(14) 


2.1 


(14) 


2.0 


(14) 


CH3OD 


1.7 


(13) 


1.7 


(13) 


2.4 


(13) 


2.4 


(13) 


N 2 


1.9 


(17) 


3.2 


(16) 


3.2 


(16) 


1.9 


(17) 
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Table 5. The molecular D/H ratios calculated in the cloud model and in the disk. For the 
disk we present two ratios (1) the midplane value and (2) N(XH)/N(XD), where XH and 
XD are the non-deuterated and deuterated forms of the molecule respectively. The ratios 
are given in the form a(—b) = a x 10~ 6 . • • • indicates that the molecular gas phase 
abundances are too low to determine meaningful D/H ratios. The midplane ratios in 
Models B and C are similar since the photodesorption included in the latter model is not 
efficient in the midplane. In Models A and D, the lack of CRH means that most species are 
frozen out in the midplane and this high depletion results in high deuteration of H3 . 



Molecule Cloud Midplane N(XH)/N(XD) 

Models A & D (Models B & C) Model A Model B Model C Model D 



Gas Species 



H 2 D+/H+ 


1.6 


(-1) 


9.9 (-1) 


9.5 (-1) 


1.9 


(-1) 


2.6 


(-1) 


8.7 


(-1) 


9.1 


(-1) 


HD+/H+ 


2.4 


(-2) 


1.6 


1.5 


2.4 


(-1) 


3.5 


(-1) 


1.2 


(-1) 


1.2 


(-1) 




^ 1 
O.l 


< 'X\ 

\- 6 > 


7& Q 


ZD. 1 






A 1 




14.6 


49.1 


CH 2 D+/CH+ 
CHD+/CH+ 


8.7 


(-2) 






1.1 


(-2) 


1.1 


(-2) 


3.8 


(-3) 


3.6 


(-3) 


3.7 


(-3) 






4.8 


(-4) 


4.7 


(-4) 


9.1 


(-5) 


8.3 


(-5) 


CD+/CH3+ 


6.8 


(-5) 






1.0 


(-5) 


1.1 


(-5) 


1.5 


(-6) 


1.5 


(-6) 


DCO+/HCO+ 


7.2 


(-2) 




12.1 


1.2 


(-1) 


9.1 


(-1) 


1.0 




1.7 


(-1) 


DCN/HCN 


3.3 


(-2) 






1.4 


(-2) 


4.9 


(-2) 


3.6 


(-2) 


3.6 


(-2) 


DNC/HNC 


1.2 


(-2) 






4.5 


(-3) 


5.9 


(-3) 


2.7 


(-2) 


2.8 


(-2) 


HDCO/H 2 CO 


4.2 


(-2) 






2.0 


(-2) 


2.3 


(-2) 


2.2 


(-2) 


2.2 


(-2) 


D 2 CO/H 2 CO 


1.0 


(-3) 






1.8 


(-3) 


2.0 


(-3) 


3.1 


(-4) 


3.1 


(-4) 


NH 2 D/NH 3 


1.3 


(-2) 






1.3 


(-2) 


6.6 


(-1) 


9.3 


(-2) 


8.6 


(-2) 


NHD 2 /NH 3 


1.6 


(-4) 






9.3 


(-4) 


5.2 


(-3) 


1.9 


(-3) 


1.6 


(-3) 


ND3/NH3 


2.7 


(-5) 






5.9 


(-4) 


1.7 


(-3) 


3.0 


(-4) 


2.3 


(-4) 


HDO/H 2 


4.4 


(-2) 






5.2 


(-2) 


2.3 


(-1) 


1.3 


(-1) 


1.3 


(-1) 


D 2 0/H 2 


2.0 


(-4) 






5.6 


(-4) 


3.3 


(-3) 


4.4 


(-3) 


4.4 


(-3) 


CH3D/CH4 


6.3 


(-2) 






1.3 


(-1) 


1.4 


(-1) 


1.2 


(-1) 


1.1 


(-1) 


CH 2 D 2 /CH 4 


3.4 


(-3) 






1.0 


(-2) 


1.2 


(-2) 


8.6 


(-3) 


6.8 


(-3) 


CHD3/CH4 


1.2 


(-4) 






3.8 


(-4) 


5.5 


(-4) 


3.8 


(-4) 


1.8 


(-4) 


CD4/CH4 


4.4 


(-7) 






4.5 


(-6) 


2.1 


(-5) 


2.0 


(-5) 


1.0 


(-6) 


N 2 D+/N 2 H+ 


7.9 


(-2) 




12.0 


5.1 


(-1) 


7.0 




5.7 




3.5 


(-1) 


CH3OD/CH3OH 


2.5 


(-2) 






1.2 


(-2) 


1.8 


(-2) 


8.0 


(-2) 


8.0 


(-2) 


CH 2 DOH/CH 3 OH 


7.3 


(-2) 






5.1 


(-2) 


5.4 


(-2) 


8.7 


(-2) 


8.7 


(-2) 


D/H 


2.1 


(-2) 


1.5 


1.0 


1.7 


(-4) 


1.4 


(-4) 


2.7 


(-4) 


2.7 


(-4) 


vlantle species 


























HDO/H 2 


2.0 


(-2) 


2.1 (-2) 


2.4 (-2) 


2.3 


(-2) 


2.7 


(-2) 


2.8 


(-2) 


2.3 


(-2) 


D 2 0/H 2 


1.5 


(-4) 


1.8 (-4) 


1.3 (-3) 


3.3 


(-4) 


1.6 


(-3) 


1.6 


(-3) 


3.4 


(-4) 


NH 2 D/NH 3 


2.4 


(-2) 


2.4 (-2) 


2.3 (-1) 


3.4 


(-2) 


2.3 


(-1) 


2.3 


(-1) 


3.4 


(-2) 


NHD 2 /NH 3 


9.2 


(-5) 


1.2 (-4) 


5.5 (-2) 


2.1 


(-3) 


5.2 


(-2) 


5.2 


(-2) 


2.1 


(-3) 


CH3D/CH4 


6.6 


(-2) 


6.7 (-2) 


2.1 (-1) 


8.8 


(-2) 


2.0 


(-1) 


2.0 


(-1) 


8.8 


(-2) 


CH 2 D 2 /CH 4 


1.9 


(-3) 


3.8 (-3) 


4.6 (-2) 


7.6 


(-3) 


4.1 


(-2) 


4.1 


(-2) 


7.5 


(-3) 


DCN/HCN 


6.1 


(-3) 


2.0 (-2) 


1.9 (-2) 


2.1 


(-2) 


1.8 


(-2) 


1.8 


(-2) 


1.9 


(-2) 



Table 5 — Continued 



Molecule Cloud Midplane N(XH)/N(XD) 

Models A & D (Models B & C) Model A Model B Model C Model D 



Table 6. A comparison of the results of Models B, C and D with the observations. The 
calculated column densities (cm -2 ) are given at t — 1 Myrs and R = 250 AU. a(b) 
represents a x 10 6 . For LkCal5 the column densities calculated from the single dish 
observations assume that the disk has a radius of 450 AU. For TW Hya the disk is assumed 
to be 165 AU in radius. For each model N(H 2 ) = 4.2 xlO 22 cm -2 . 



Molecule 


Model B 


Model C 


Model D 


DM Tau 1 


LkCal5 


TW Hya 3 














Interferometer 2 


Single dish 3 




CO 


1.2 (17) 


2.9 


(17) 


2.7 (17) 


5.7 (16) 


1.7 (18) 7 


1.9 (16) 


3.2 (16) 


HCN 


1.8 (10) 


3.3 


(13) 


3.3 (13) 


2.1 (12) 


2.4 (13) 


1.8 (12) 


8.5 (12) 6 


HNC 


3.7 (10) 


3.2 


(12) 


3.1 (12) 


9.1 (11) 


< 5.4 (12) 




< 1.4 (12) 


CN 


3.8 (11) 


3.6 


(13) 


3.6 (13) 


9.5 - 12 (12) 


9.7 - 25 (13) 


1.5 (13) 


6.6 (13) 


CH 3 OH 


6.5 (5) 


1.5 


(11) 


1.5 (11) 




7.3 - 18 (14) 


< 7.1 (13) 


< 1.1 (13) 


H 2 CO 


3.3 (11) 


3.9 


(12) 


3.9 (12) 


7.6 - 19 (11) 


7.2 - 19 (12) 4 


7.1 - 51 (11) 


< 8.0 (11) 


HCO+ 


3.3 (12) 


3.2 


(12) 


2.8 (12) 


4.6 - 28 (11) 


1.5 (13) 


3.3 (11) 


8.5 (12) 6 


C 2 H 


3.6 (12) 


7.2 


(12) 


7.2 (12) 


4.2 (13) 








N 2 H+ 


2.3 (10) 


2.8 


(10) 


1.1 (10) 


< 7.6 (11) 


3.1 (13) 7 


< 1.4 (12) 


< 1.0 (13) 
3.0 (ll) 6 


DCO+ 


3.0 (12) 


3.2 


(12) 


4.8 (11) 






< 2.9 (11) 


DCN 


8.9 (8) 


1.2 


(12) 


1.2 (12) 






1 (13) 8 


< 4.0 (10) 


HDO 


9.3 (9) 


2.0 


(14) 


2.0 (14) 


1.6 (13) 




2-7 (14) 8 




H 2 D+ 


8.2 (11) 


7.7 


(11) 


4.8 (11) 


8.8 (12) 5 




< 2.0 (14) 5 


1.3 (14) 



Reference s. — (1) lAikawa et all (120021) (derived fromlDutrev et all (Il997h). (2)lQil (l200ll) (3) iThi et al 
(2004). (4) lAikawa et all |2003l ) ~5) ICeccarelli et all |2005l ). (6) Ivan Dishoeck et all |200a ). (7) foi et al. 
(|2003l ). (8) bil (j200ll )~ 
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Fig. 1. — The density and temp erature distributions used in this paper. These are from the 
models of Id'Alessio etaD J200lh . with M = 10~ 8 M Q yr" 1 , M* = 0.7 M , T* = 4000 K, R* 
= 2.5 R Q . The grain size distribution is given by n(a) oc a~ 3 5 for 0.005 < a < 0.25 /zm, 
where a is the grain radius. The surface density So = 25 gem" 2 at 10 AU and varies as 1/R 
for R > 10 AU. The mass of the disk is 0.063 M out to R = 400 AU. 
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Fig. 2. — Schematic of the chemical structure of a disk, showing how it can be divided into 
three layers. In the midplane the high densities and low temperatures ensure that most 
molecules are frozen out onto the surfaces of the dust grains. At the surface, the high UV 
field dissociates molecules into their constituent atoms and ions. Between these two lies a 
molecular layer where the observed emission originates. 
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Fig. 3. — The fractional abundances calculated at R = 250 AU at a time of 1 Myrs for Model 
B as a function of z. 
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Fig. 4. — The fractional abundances calculated at R = 250 AU at a time of 1 Myrs for Model 
C as a function of z. 
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Fig. 5. — The fractional abundances as a function of height, z above the midplane for Models 
B and C. In Model B, the sharp decrease in abundance of CO and CH 4 at z ~ 60 - 70 AU, 
and of H2CO and isotopomers at z ~ 120 AU is due to the onset of thermal desorption. 
The carbon atoms produced in the gas phase when CO is released and then broken up are 
reaccreted on to the grains between 60 and 140 AU where they form hydrocarbon ices. At 
z > 140 AU, most accreted carbon atoms are desorbed before they can react and hence the 
abundance of hydrocarbons in the ice falls off. In Model C photodesorption ensures that all 
molecules are removed by from the grains by z = 80 - 100 AU. 
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Fig. 6. — The radial column density distributions calculated in Model B. 
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Fig. 7. — The radial column density distributions calculated in Model C. 
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Fig. 8. — The D/H ratios calculated from the ratio of the column density of the deuterated 
species to the column density of the non-deuterated species, as a function of R. Model B is 
shown in the solid line, with Model C as the dashed line. 
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